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1. MATERIALS ENGINEERING 


Senior Members 
Dr. Manfred Altman 
Dr. Solomon Pollack 



1 . 1 THERMAL DIFFUSIV1TY OF MIXTURES 


Senior Investigators: Dr. M. Altman and Dr. K. Sreenivasan 
Graduate Student: S.J.Amir 


OBJECTIVES 


To measure the thermal diffusivity of two types of mixtures: 

1. solid mixtures in which the components of the 
mixtures are both in solid state 

2. solid-gas mixtures in which the solid has a certain 
amount of porosity containing gases like helium or 


PREVIOUS ACCOMPLISHMENTS 

The thermal diffusivities of the following components and mixtures 
were reported: 


1: 

CaF 2 

INDEC-SR and SR 

-9 

2. 

BaF 2 

INDEC-SR-9 


3. 

MgF 2 

INDEC-SR-13 


4. 

CaF 2 + BaF 2 

(1 to 1 weight) 

INDEC-SR-11 

5. 

CaF 2 + BaF 2 

{1 to 2 weight) 

INDEC-SR- 11 

.6. 

Ca F 2 + Ba F 2 

(2 to 1 weight) 

INDEC-SR-11 

7. 

CaF 2 +MgF 2 

(1 to 1 weight) 

INDEC-SR-12 , SR-13 

8. 

CaF 2 + MgF 2 

(1 to 2 weight) 

INDEC-SR-13 


The thermal diffusivity was measured over the range of 100 
1100 °c using a vacuum furnace. 


o 


c to 
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PROGRESS IN PAST PERIOD 


The success of research on determination of thermal diffusivity of 
porous media at high temperatures depends on the preparation of sample 
with high porosity, high purity, reasonably uniform pore diameter and 
uniform density. In the past period many methods of sample fabrication 
were tried, and a method to make a sample with definite porosity- and 
other characteristics was developed and finalized. Details are given - 
elsewhere in this report. 
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1.2 TRANSIENT TEMPERATURE RESPONSE OF CYLINDRICAL, 
COMPOSITE ENERGY STORAGE DEVICES 


Dr. M. Altman, Director 
Dr. K. Sreenivasan, Post Doctoral Fellow 


OBTECTIVES 

To develop the analytical solution for the temperature distribution 
in a cylindrical, composite energy storage device. 


PREVIOUS ACCOMPLISHMENTS 
None 

PROGRESS IN PAST PERIOD 

The transient solution to the heat conduction equation in a composite 
cylinder was developed assuming perfect thermal contact at the interfaces. 

The solution was used to obtain the temperature response of the energy storage 
device. Details are shown on pages A-i-6 through A 1- 11 . 
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1.3 EXPERIMENTAL DETERMINATION OF THE THERMOELECTRIC 
PROPERTIES OF GRAPHITE ALLOYS 

Senior Investigator: Dr. S. R. Pollack 
Graduate Student: J. J. Curry 

OBJECTIVES 

To determine the usefulness of quaternary tetrahedral compounds 
as thermoelectric generator materials. 

PREVIOUS ACCOMPLISHMENTS 

Reproducible powder metallurgical techniques were developed for 
fabricating these materials. Electronic equipment was built for measur- 
ing thermoelectric power, electrical conductivity, and Hall coefficient as 
a function of temperature. These measurements were made, and although 
a high thermoelectric power is obtained, the low value of electrical con- 
ductivity prevented usefulness as device material. 

PROGRESS IN PAST PERIOD 

Analysis of data continues seeking an understanding of the low 
conductivity.. It was hoped that it could be improved if the mechanism 
were understood. Unfortunately, we believe the mechanism is alloy 
scattering, and in order to circumvent this problem, materials fabrication 
techniques superior to all those tried during the course of this work will 
have to be developed. 

This work was presented at an international conference, and is in 
preparation for publication. 
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1.4 SUPERCONDUCTIVITY IN EVAPORATED TUNGSTEN FILMS 


Senior Investigator: Dr. S. R. Pollack 
Graduate Student: S. Basavaiah 


OBTECTIVES 

To determine the origin of the Tc enhancement in evaporated 
tungsten films, and to study the properties of the superconducting films. 


PREVIOUS ACCOMPLISHMENTS 

The origin of the Tc enhancement was traced to the presence of 
(3 -tungsten (A15 structure) which was stabilized during deposition by the 
presence of oxygen . Control of Tc can be achieved by controlling the 
oxygen partial pressure, thereby controlling the relative amounts of 0 to 
a tungsten. Using tunneling techniques, the energy gap of 8 -W was 
measured as a function of Temperature, and found to be adequately de- 
scribed by B . C . S . Theory , 


PROGRESS IN PAST PERIOD 


This work was completed and two papers were published. The 
work is continuing with other materials under sponsorship of ARPA. 
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1.5 BOILING THERMIONIC CONVERTER 
Senior Investigator: Dr. Alan M. Whitman 


OBTECTIVES 

To study the feasibility of a liquid metal collector thermionic con- 
verter stabilized by gravity. 


PROGRESS IN PAST PERIOD 

An experiment was performed with a low boiling point, non-conducting 
fluid, wherein a stable vapor-liquid interface was maintained with arbitrarily 
small vapor film thickness. In addition, a thermodynamic analysis of the 
system was performed which yielded a set of equations suitable for the design 
of such a device . 
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2. PLASMA ENGINEERING 


Senior Members 
Dr. Louis Girifalco 
Dr. George L. Schrenk 
Dr. Hsuan Yeh 



2.1 MHD FLOW IN AN ANNULAR GAP 


Senior Investigator; Dr. I. M. Cohen 
Graduate Students: W. C. Frazier, J. P. Hanson 

QBTECTIVES 

To study the fluid mechanics of rotating MHD flows . 

PREVIOUS ACCOMPLISHMENTS 

* 

Initial experiments were carried out by T. Ebtekar on an old electro- 
magnet found in the laboratory. This established the orders of magnitude 
and trends to be expected in the primary variables . With this information, 
a new magnet was ordered and more careful experiments were prepared. 

PROGRESS IN PAST PERIOD 

** Measurements of the tangential component of velocity as a function 
of the current driving the magnet and the current through the device were made 
using a pitot tube. These measurements were taken at a point halfway between 
the inner and outer electrode, and 0.95 in. below the top electrode. In ad- 
dition, two velocity profiles from the top electrode to the point 0.95 in. below 
it were taken at constant device and magnet currents . These profiles showed 
the tangential component of velocity to be independent of the axial position 
over this distance (about one sixth of the total length of the device) . 

The results of the measurements taken as a function of magnet and de.- 
vice current are shown in the attached figure. The pressure differences shown 
are between the total pressure at the pitot tube, and the static pressure at 

*T . Ebtekar, "Flow of a Conducting Liquid in an Annular Gap", Ph.D. disserta- 
tion, University of Pennsylvania , May, 1968. 

** John P. Hanson 
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the upper insulator at the same radial position as the pitot tube, corrected 
for the difference in gravity head. The velocity calculation indicator on 
the figure neglects the magnetic pressure, assumes the r and z components 
of velocity are small, and that the axial static pressure gradient is only 
that due to gravity. 

The curves in the figure show that the velocity goes as the square 
root of the device current, if the assumptions concerning the pressure 
gradient at secondary velocity are born out by further experiments and/or 
analysis. The line marked "Critical Line" is an overlay of data obtained 
in other experiments, and are explained more fully in Mr. Frazier’s section 
of the report immediately below. 

To determine the nature of the MHD flow, i.e. turbulent or laminar, 
several devices were constructed and placed in the flow generated in a 
cylindrical gap. The most promising is a streamlined shape that approximated 
within 10% NACA 23024 made from stainless steel. Span is 0.50" , chord 

0. 25" . A torsion arm was inserted at 0.065" from the leading edge. Angle 
of attack changes (flutter) was then carried by this torque arm to a light 
interruption system that measured the changes in angle of attack by propor- 
tional voltage from a photocell. Since the major part of the moving mass is 
in the NACA shape and bearing drag is very low, the system is capable of 
following fluid incident angle changes on the order of 30 to 50 cycles per 
second. Computations from the lift -drag data associated with this NACA 
shape indicate a natural flutter at 40 cps . 

Inserted in the flow at 1/3 the distance from the top of the device, 
and at approximately 2/3 the distance from the inner electrode, the shape 
was found to flutter at higher frequencies as the tangential velocity was 
increased. At low velocities, the device was stable and did not flutter — 
see included plot. 

Additional experiments have shown thereto be two vibration modes; 

1. e. , the shape flutters about a zero angle of attack to r tangential and at 
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an angle of about 30° to this zero line. Further modification of the device 
will be accomplished to measure these modes accurately. 


SeS Cell 
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2.2 BASIC SURFACE INVESTIGATIONS 


G . L . Schrenk , S . P . Sharma 
OBTECTIVES 

To study the influence of high fields on surface charge distributions . 


PREVIOUS ACCOMPLISHMENTS 

Research is in progress to understand the structure of field ion 
micrographs . In particular, we are trying to understand the varying average 
intensities of different regions and also the varying intensity of various dots. 
The regions of varying intensity are rather sharply defined, and are more 
pronounced in some metals than in others; e.g. , in platinum they form one 
of the most prominent features of the micrograph, while in tungsten, they 
are barely discernible. 

Platinum, as imaged by helium, has been extensively studied. 

Definite positive effects from the shape of the platinum Fermi surface have 
been obtained; these results can explain a large part of the structure of the 
platinum micrograph. These results are limited by the lack of knowledge of 
the Fermi surface of platinum. This calculation is the first known calculation 
of high field tunneling probabilities that takes into account the detailed 
structure of the Fermi Surface of the metal and was the Ph.D. dissertation 
for Steve Fonash. This calculation was summarized in SR-14. 

PROGRESS IN PAST PERIOD 

Research is currently underway to understand the varying intensity 
of the individual dots on a field ion micrograph. We are studying such 
basic questions as which surface atoms are imaged in a field ion microscope 
and why . 

A quantum mechanical model of the tunneling process has been 
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constructed; this model treats the variation of surface potential in the plane 
of the surface as well as the variation perpendicular to the surface. Exten- 
sive calculations are in progress. Preliminary results indicate that the 
amplitude of the variation of the .surface potential in the plane of the surface 
can have a very significant effect on the tunneling probability. These 
results are the basis of the Ph.D. dissertation of S . P. Sharma . 
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2.3 A MATERIALS STUDY OF SILICON SOLAR CELLS 


Senior Investigator: S. R. Pollack 
Graduate Student: W. H. Becker 


OBJECTIVES 

To study the structure and properties of contacts to Si with special 
attention given to titanium -silver contacts used in solar cells in conjunc- 
tion with the space program . 

PREVIOUS ACCOMPLISHMENTS 

A survey of the field was made in order to determine the state of 
understanding. Unfortunately not much work has been done regarding the 
materials science aspects of the Ti-Ag contacts. 

PROGRESS IN PAST PERIOD 

Production solar cells from several manufacturers were obtained. 
Structure and Phase identification work has been started on "good" cells 
and failures usir?g X-ray diffraction and electron microscopy . Techniques 
were developed for preparing electron microscopy samples of material under 
the -contact blisters that occur during failure . Electron diffraction data and 
transmission electron micrographs have been taken and are currently being 
analyzed. 
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2.4 THERMOELECTRIC PROPERTIES OF GRAPHITE COMPOUNDS 


Principal Investigator: L. A. Girifalco 
Post-doctoral Fellow: S. Sachidanandam 

Graduate Student: T . Montelbano 

OBJECTIVES 

The objective of this work is to examine lamellar alloys of graphite 
in order to understand the relation of their structure to thermoelectric prop- 
erties with a view to the possible use of such compounds for thermoelectric 
applications . 

PREVIOUS ACCOMPLISHMENTS 

An extend Kronig -Penney model was used for calculations of the 
electronic band structures of lamellar compounds. The results showed how 
parameters such as crystal potential and lattice parameters affect the 
structure, and the semi -conducting or metallic properties of such compounds. • 

Theoretical work on acoustic phonons in lamellar compounds was 
performed, and methods of solving the lattice dynamics problem were developed. 

Specimens of barium-graphite alloy were prepared. This is a new 
compound, and its preparation involved many experimental difficulties. The 
most concentrated compound C 0 Ba, is hard and brittle, and has a gold 

* O 

metallic color . 

PROGRESS IN PAST PERIOD 

The lattice dynamics work was completed and applied to an analysis 
of the phonon conductivity of graphite alloys . The phonon conductivity was 
found to be reduced significantly below that of pyrolitic graphite. However, 
the increases in the electrical conductivity and thermoelectric power, as 
reported from measurements by other- workers on metal -graphite lamellar 
alloys, give a thermoelectric figure of merit two orders of magnitude below 
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the commercially viable values . Considering the wide range of values the 
various samples of graphite are found to have for the electrical parameters, 
there is reason to believe that purer samples could have much higher figures 
of merit . The electrical properties are being studied with a view to explor- 
ing this possibility. 

A paper entitled "Acoustic Phonons in Some Metal -Graphite Lamellar 
Compounds" is in press with the Physical Review. 

The experimenta 1 work on barium graphite is continuing. 
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Senior Members 


Dr. John O'M. Bockris 
Dr. Leonard Nanis 





3.1 OVERPQTENTIAL TRANSIENTS ON THE ROTATING DISC 


Senior Investigator: Dr. E. Nanis 
Graduate Student: Irving Klein 

OBJECTIVES 

The present problem is an attempt to treat in a rational way a typical 
electrode for which the mass transport is controlled to permit prediction of 
electrochemical effects. The rotating disc geometry was found to be ideal 
for this purpose, and may be used for actual battery electrodes and electrolytes. 
The present study will accomplish a long desired treatment of combined ac- 
tivation and concentration overpotential transient behavior. Engineering 
extension to fuel cells and batteries may be obtained by knowledge of the 
local Nusselt numbers which control mass transport at working electrodes. 

PREVIOUS ACCOMPLISHMENTS 

A mathematical model has been constructed for the transient redox 
overpotential decaying from steady state following interruption of galvano- 
static conditions. The model includes initial conditions for the individual 
concentration profiles of an oxidized and reduced species and the time 
dependent Navier-Stokes equation as applied to a rotating disc geometry. 

Also, a time dependent functional expression relating overpotential of the 
total redox system to the concentration of the individual component species 
has been, derived. A rotating disc system has been designed for experimental 
verification of the analytic treatment, and was made operational. The ap- 
paratus was used to examine the steady state behavior of potassium ferri- 
cyanide and potassium ferrocyanide redox couples; i.e. , the variation of 
limiting current densities with speed of rotation, and confirmation of 
theoretical hypotheses was achieved. 
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PROGRESS IN THE PAST PERIOD 


Transient overpotential decay data for the potassium ferricyanide- 
potassium ferrocyanide redox couple have been obtained. Comparison 
of the data with theoretical values indicates a good correlation and serves 
to verify the mathematical analysis upon which the present research is 
based. Transient overpotential data have been obtained as a function of 
the following experimental parameters; bulk concentration of the reacting 
species, bulk concentration of the supporting electrolyte, the rotational 
•speed of the disc, fraction of the limiting current prior to current shutoff. 
Deviations from the theory occur only when the bulk concentrations of the 
reacting species are less than 5xl0 4 gmole/liter. The observed anomalies 
can be explained by considering the possible complexing of the ferricyanide 
ions with the large excesses of potassium ions created by the presence of 
the supporting electrolyte, potassium hydroxide. 

Work involving the study of transient overpotential buildups at 
rotating disc electrodes has been initiated. A mathematical model analogous 
to the one developed for transient overpotential decays has been derived. 
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3.2 CURRENT AND POTENTIAL DISTRIBUTION IN CYLINDRICAL 

GEOMETRIES: ENGINEERING APPLICATION TO FUEL CELL DESIGN 


Senior Investigator: Dr. L. Nanis 
Graduate Student: Wallace Kesselman 

OBTECTIVES 

The purpose of this research is to derive simple mathematical ex- 
pressions useful for engineering evaluation of current and potential distri- 
butions for common geometries, as encountered in fuel cells or battery 
design. 


PREVIOUS ACCOMPLISHMENTS 


Current and potential distributions on a disc under conditions of 
activation polarization were worked out and described in terms of the 

STL 


throwing power parameter, k/a 


5 J 


Similarly, the distributions 


were solved for the cylinder-versus -plane electrode configuration for the 
primary distribution. 


PROGRESS IN PAST PERIOD 

The cylinder-versus -plane configuration was further analyzed to in- 
clude the effect of activation polarization on the cylinder. Current distribu- 
tions were obtained in agreement with the empirical predictions for the 
system with overpotential effects considered. 

Furthermore, experimental measurements were taken on this latter 
configuration to verify the distributions obtained in the last period. Good 
agreement between theory and experiment was found for the primary current 
distribution . 
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3.3 FOAMING ELECTROLYTE FUEL CELL 


Senior Investigator: Dr. L. Nanis 
Student: Frank R. McLarnon 

OBJECTIVES 

The diffusion of reacting species provides an ultimate limitation in 
all electrochemical processes. For fuel cell applications, dissolved electro- 
active gases may be rapidly dissolved in a very thin film of electrolyte. 

Such a structure is provided by a foamed electrolyte and may be useful in 
practice if the rate of transport is sufficient and the circuit resistance is 
low. Verification of this possibility is the goal of present laboratory testing 
of cell geometries, surfactants' and electrolytes. 

PREVIOUS ACCOMPLISHMENTS 

It has been demonstrated that electrode reaction can be supported 
through a foamed electrolyte. 

Work was completed on the design, construction, and testing of a 
cell which can be readily disassembled and cleaned, and which minimizes 
migration of oxygen into the anode compartment . 

PROGRESS IN PAST PERIOD 

Several foam producing compounds have been studied with hydrogen 

in 4N H„SO . oxidized on platinum (smooth) anodes. Current densities of 
^ ^ 2 • 

up to 15 mA/cm have been achieved. Continuous operation for 4 hours 

resulted in an increase of anode overpotential from 24 to 34 millivolts at 

2 

a current density of 0.9 mA/cm . 
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3.4 HIGH PERFORMANCE IODINE CATHODE 

Senior Investigators: 

A. P. Saunders Dr. H. Wroblowa 

OBTECTIVES 

To investigate and develop the porous cathode previously reported 
under "Primary and Secondary Batteries using Charge Transfer Complexes" . 


PREVIOUS ACCOMPLISHMENTS 

A means was found of obtaining high power densities and high 
faradaic efficiencies for fuel cell consuming Mg and I ^ , by means of a 
porous cathode; the I being in a concentrated solution and flowed under a 

c* 

pressure gradient through the finely-pored graphite cathode. 

PROGRESS IN PAST PERIOD 

A current potential relationship has been studied as a function of 
the pressure gradient across the cathode/ The faradaic efficiency at the 
maximum steady state current at a given'pressure is 100% (within the ex- 
perimental error of± 5%). The power loss at 95% efficiency of Ig consump- 
tion at the cathode is low, the loss increasing with further increase in 

faradaic efficiency. The design difficulties are discussed. The optimum 

-2 

power density of a Mg I^ fuel cell is evaluated at 0.3 watts cm , the 
limiting electrode being the metal anode. 
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THERMAL DIFFUSIVITY OF MIXTURES 


Senior Investigator: Dr. M. Altman and Dr. K. Sreenivasan 
Graduate Student: S. J. Amir 

THEORY AND GENERAL METHOD OF APPROACH 


In principle the thermal diffusivity of solids can be measured in any 
apparatus that provides the requisite boundary conditions to the equation 
for heat conduction in an isotropic material without heat sources and for 
constant physical properties. 


u _ 

7>t 


oi V 2 T 


where Ofis the thermal diffusivity defined by 

The solution of the equation for transient conduction in a simple geometry 
will give the thermal diffusivity from the measured temperature distribution 
as a function of time. The thermal conductivity can then be calculated from 
a knowledge of the specific^heat and density of the solid. Such a method is 
usually referred to as the "transient" method. It 'should be noted that no heat 
flux measurements are necessary. 

Steady state methods rely on the solution to the equation 



V a T i O 


together with information on the heat flux and temperature at the boundaries. 

The thermal conductivity is directly obtained in this method. The thermal 
diffusivity can be calculated from a knowledge of the specific heat and density. 
Such a method is usually referred to as the "steady state" method. 

Accurate determination of heat flux at high temperature is extremely 
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difficult. Hence, the quasi-method was chosen for the measurement 
of thermal diffusivity. 

The method can be described as follows: Fundamentally the 

difference between the temperature at the surface and at the center of a 
cylindrical specimen is measured for a constant rate of surface temperature 
rise. The thermal diffusivity is then calculated using the diameter of the 
specimen, the rate of rise of the surface temperature and the temperature 
difference. The boundary condition on which this method is based can be 
easily realized experimentally. No heat flux need be measured. This is a 
tremendous advantage as accurate measurement of heat flux at high temper- 
atures is very difficult. The mathematical details of this method are given 
in references 1 and 2 . 

Thermal Diffusivity of Mixture Samples 
In the void fraction of two constituent species A and B are and v g 
the thermal diffusivity of the composite sample can be predicted only if 
its thermal conductivity and the specific heat per unit volume are known. 
The composite specific heat per unit volume can be obtained by 

(f<?\ = T'a (? c r) A 4 v * It >c f)t 

f - Density - lbs /c.ft V A void fraction of A 

Cp Specific heat - /jrw J lb a f V voi< ^ fraction of B 

Subscript T - composite 
A - Species A 
B - Species B 

In the case of thermal conductivity of the composite one can only 

i 

calculate the upper and lower bounds. The two bounds are given by 

/<r, = ^ k A + z/g kg 

kr z * k& - 

v k f<4 


(2) 

(3) 
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The bounds can be improved from a knowledge of statistical' 
correlation functions of the mixture samples. The aim of this proposed 
investigation is to determine the correlation functions experimentally and 
study how to improve the above mentioned bounds 

Thermal Diffusivity of Porous Samples 

In this case one of the constituents is gas and its properties have to 
be used in the equations (1), (2) and (3). If the material can be considered 
structurally as a solid matrix with isolated pores the thermal conductivity 
is given to a first order approximation as . 

jc v + Up + 4 E o- dp T 3 J) 

where subscripts s, p and g refer to the solid, pore and gas', respectively 
and 

E - Total emissivity of the inside surface of the pore 

T - Absolute temperature 

d - width of the cavity 
P 

The above equation is obtained by considering that the gas in pore 
does not participate in the radiation transfer of heat. It is the purpose of 
this investigation to recast the above formula assuming that the gas participates 
in the radiation heat transfer^ 

Making the Specimen 

In the previous investigation on solid specimens no difficulty was 
experienced in obtaining them. The Lucalox samples were bought, and the LiF 
samples were made by melting LiF sample in a carbon mold heated by high 
frequency eddy currents in an induction furnace. However, it is difficult to buy 
samples with a definite porosity or a mixture sample of different ratio of the 
constituent species. 

EXPERIMENT: 

The experimental program consisted of preparing samples of aluminum 
oxide, magnesium oxide, titanium oxide and zirconium oxide. The main 
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interest was for aluminum oxide since its conductivity and emissivity are 
well established. Titanium oxide, besides aluminum oxide was also found 
amenable to specimen fabrication. Magnesium oxide was, however, found 
very hard to work with while there was partial success with zirconium oxide. 

All samples were prepared by incorporating Napthelene powder in the 
oxide powder. Just before use Napthelene powder was passed through a 
screen with opening of 0.42 MM. This would induce uniformity of the pore 
diameter. The mixture of oxide and Napthelene was then mixed with organic 
liquid binder. The liquid binders were metal resinates with metal content up 
to 10% . 

The aggregate of binder, Napthelene and Metal oxide was then pressed 
in a die at a pressure of 3200 PSI. For Aluminum and Zirconium the die was 2" 
diameter while for Titanium and Magnesium it was 2 ~ " . The bigger diameter 

O 

size used was to account for the shrinkage taking place in those oxides. 

After pressing, the samples were dried for an extended period of time at 70-80 c 
until no Napthelene odor could be noticed. The samples were then fired at 
1280 °c in air for about 24 hours so that the organic material carbonizes and the 
metal from the resinate should form its respective oxide. 

The pure oxide sample now is fired in vacuum at temperatures 80% of the 
melting point of the metal oxides for four or five hours in standard manner. 

After firing, polished sections were examined which showed that pores were 

w 1 " 
roughly sperical or ellipsoidal in shape. Samples (2" dia. and 1— deep) with 

up to 65% porosity were successfully made. 
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TRANSIENT TEMPERATURE RESPONSE OF CYLINDRICAL, COMPOSITE 

ENERGY STORAGE DEVICES 


M. Altman, Director 
K. Sreenivasan, Post Doctoral Fellow 


This paper describes an analysis of the transient temperature response 
of cylindrical, composite energy storage devices whose thermal conductivity, 
density and specific heat are independent of temperature. Solutions are 
presented for a three-layer composite cylinder with heat fluxes specified at 
the inner and the outer faces. _The conventional method of separation of 
variables is used. It is shown that the customary weighting function g for the 
Bessel functions, has to be modified to g, n respectively for the 

three layers to make the eigenfunctions orthogonal in the region considered. 
Expressions are derived for the amount of energy available for storage and 
efficiency of storage. 
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The geometry of the system considered is described in Fig. 1. The 
composite cylinder receives heat from outside and this heat is conducted inwards 
to the core where it is utilized. In this process it heats up the intermediate 
storage material. The storage material has to be properly sheathed and increase 
in its enthalpy during any interval is a measure of the energy available for 
storage. 

Conceptually, the physical analogue would correspond to a situation where 
the regions (1) and (3) would be sheathing material containing the energy storage 
material in region (2). The inner core will contain the energy conversion devices 
and the outer surface of region (3) will receive heat ‘from outside. 


The aim of this analysis is to determine the following: 

(1) The optimun geometry for the device 

(2) The influence of the. thermal conductivity, thermal diffusivity specific 
heat and density of the various layers , on the over all performance 

(3) Variation of temperature and heat flux at every point in the device. 

To demonstrate the method of analysis it was decided to consider a four- 
layered composite cylinder. The following assumptions were made. 

(1) Heat flow in the radial direction only. 

(2) Perfect thermal contact at the interfaces. 

(3) Physical properties do not vary with temperature. 

(4) No internal heat generation. 

(5) The material in the four layers is isotropic and homogeneous. 

The boundary conditions considered in this analysis are for the case where 
the heat fluxes at the inner and outer boundaries of the composite cylinder are 
specified. 

This is usually the case in practice where the heat fluxes are specified and 
one wants to predict the corresponding variation. The following dimensionless 
variables are employed. 
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where i=l, 2 and 3 for the three regions. In terms of these variables the 


governing equations are 

= i' v *t 

7>e c* 3 '• 

^ - ^2 v'T 

d© 7* v 

= V X T, 
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•i ^ 
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with the initial condition at 0= 0, %. C*j, o) = T 3 fy,op=. O 

and the boundary conditions » _ ft . ^5 - 

jw 7_ *• : * k - 

The matching conditions at the interface are given by 

k, h = 1 I* 

i-*; = r s f5,*; 

it. ^ 7i 5 ^3 

3*; 37 

Unavailability of space prevents us from presenting the details of the solution. 
They are found in Reference 1. The final solution can be written for the three 
regions as 

TCI.*) = + s: £) £,£.-) <- C, + 

T.c.ej = A(e + ^ £ ,j 

A fe + -+ +c 3 t Xi. (ij 

z ^fi- Hi A) 

feiTis V <>• J 
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and ^ s are zeros of the secular determinant in the equation 
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The temperature variation at the inner surface is obtained by substituting r»= — 

d 

in equation (5a). The heat flux is given by 
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where the temperature derivatives can be obtained by differentiating equations (5), 

Another parameter of interest is the total heat energy available for storage. This is 

C t) 

given by the difference in heat flux crossing the interfaces r\ - — and ti 
* 




where the derivatives are calculated as indicated earlier. As F* represents a rate 
let us now obtain the total energy. If the system is in the bright period for a time 
t c the energy per unit area impinging on the outside Is 




= 

P> 3. 


Total energy available for storage is given by 

E-. = %*Jt = £ e Jr a 
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To assess the performance of the storage device it is convenient to define a 
storage efficiency by the ratio, F&wi / E ■ or 

, e . . ^ 

yi = - f f f 1^7 - f S 7 , ? 

This analysis clearly demonstrates how all the relavant parameters of 
thermal storage device can be calculated. The great advantage of this method 
of analysis lies in the fact that no horrendous inversions are involved and that 
all the calculations, though painstaking, are possible. 
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THERMOELECTRIC PROPERTIES OF GRAPHITE COMPOUNDS 


Principal Investigator: L. A. Girifalco 
Post-doctoral Fellow: S. Sachidanandam 

Graduate Student: T . Montelbano 

PREPARATION 

The barium -graphite compound C Ba has been synthesized by dif- 

O 

fusion anneal. The optimum condition for preparation of this new compound 
is heating barium and pyrolytic graphite in vacuum at 575 C for 400 hours. 
At higher temperatures the formation of barium carbide competes for avail- 
able barium and yields an inhomogeneous composite. At a lower tempera- 
ture, 500°C for example, there is approximately 50% less compound forma- 
tion after 1600 hours of anneal. 


STRUCTURE 

The lattice parameters of C 0 Ba have been determined by Debye- 

O 

Sherrer techniques . The use of thin walled capillaries allowed for all X-ray 
work to be performed in vacuum. In the " a " -direction the parameter of C^Ba 

o 

is 2.49A while that measured for the pyrolytic graphite used is 2.45A. The 
distance between carbon layers in the “c" -direction is 5.28Afor compound 
and 3 . 35A for graphite . 

Macroscopically the expansion along the c-axis is manifest by a 
doubling of height. No apparent change is noted in the a -direction. It is 
important to note that mechanical integrity of the compound is maintained 
upon formation . 
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ELECTRONIC PROPERTIES 


Resistivity measurements have been performed in vacuum in the 
temperature range from room temperature to 600°C. Preliminary results 
are shown in Figure 1 . 


to 



Figure 1. Resistivity of C g Ba in temperature range 
from room temperature to 600*0. 

The transition from semi-conducting to metallic type behavior makes this 
compound very interesting for electronic study. Changes in resistivity are 
attributed to a loss of barium from the C Ba . This would involve the forma- 

• O 

tion of less concentrated compounds and concurrent structure alterations. 
Work is continuing to elucidate these changes. 

Investigation of thermoelectric properties in the same temperature 
range is planned. 
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CURRENT AND POTENTIAL DISTRIBUTION IN CYLINDRICAL 
GEOMETRIES: ENGINEERING APPLICATION TO FUEL CELL DESIGN 

Senior Investigator: Dr. L. Nani's 
Graduate Student: Wallace Kesselman 


By means of the bilinear transformation of coordinator, the cylinder- 
plane configuration can, by suitable choice of constants of the transforma- 
tion, be reduced to one of concentric circles for which, as demonstrated last 
quarter, the current distributions can be solved and expressed in the original 
geometry in terms of simple analytical expressions. This method has been 
extended to the case where overpotential exists. In this latter case, after 
transformation of coordinates , the LaPlace equation must be solved by sep- 
aration of variables, and the results obtained in this geometry must then be 
made to correspond to the cylinder-wall geometry. The results obtained can 
be compactly plotted as front -to-back -face ratio of current densities, and 
are shown in Fig. 1 as a function of throwing power parameter. 

The expression for the primary current distribution on the cylinder 
in this geometry 


J/Jave = 


2 , 2 . 

P /a - 1 

p/a + sin -0- 


where p = distance from wall to center of cylinder 

a = radius of cylinder 

■6-= angular position on cylinder as measured from 
the horizontal in a counter-clockwise direction 

was verified by means of an experiment in which one wall of a pan with 
vertical walls was made as electrode whose counter-electrode was a glass 
pipe with nickel foil wrapped around it, leaving a gap for a segment in- 
sulated from the main electrode. Current was then measured through this 
segment. Agreement was achieved as shown in Figures 2 and 3 . 
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FOAMING ELECTROLYTE FUEL CELL 


Senior Investigator: Dr. L. Nanis 
Student: Frank R. McLamon 

The improved design foam fuel cell (see Sec. 3.3 INDEC SR-14, 

June, 1968) has been used for anodic oxidation studies of hydrogen in 
4N H^SO^ on bright, smooth platinum electrodes . The addition of 0. 1 
weight percent GAFAC RE-610 to the electrolyte sustained currents larger 
than the 0.6 mA/cm associated with vigorous stirring by hydrogen bubbles. 
This was taken as evidence of participation of the RE-610 additive in the re- 
action for totally immersed electrodes . The foam produced in this way per- 
mitted a higher current density (3 mA/cm ) to be reached before excessive 
polarization due to diffusion limitation. Studies of GAFAC RE -6 10 were 
discontinued because the unwanted side reactions involving the additive 
masked true foam structure effects on the limiting current density. The 
experience with GAFAC RE-610 prompted studies of background currents in 
foams produced by various additives for stirring by helium. In addition, 
electrolyte with and without additives was studied for behavior of the 

current -overpotential relation for several stirring conditions produced by 

\_ 

regulating the hydrogen flow rate. Equal currents for solution with and 
without additive (for the same stirring condition) were taken as an indica- 
tion of non-participation by the additive. As a futher check, current- 
overpotential curves for the liquid phase stirred by helium alone clearly 
demonstrated the absence of hydrogen reaction, since the region of oxygen 
overpotential was reached when the current density exceeded 0.05 mA/cm 

i.e. , much less than the liquid phase hydrogen limiting current density of 
2 

0.6 mA/cm for comparable stirring . With the above background tests, as 
a guide, electrolyte foamed by the addition of 0.1 wt. percent Textilana 
FM 1001 permitted the attainment of a limiting current of 15 mA/cm . 
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Preliminary life tests with this additive have shown a 2 mV/hr increase of 

2 

anode overpotential after 4 hours of operation at 0.-9 mA/cm , chosen to 
correspond to an initial polarization of 24 mV. 
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HIGH PERFORMANCE IODINE CATHODE 


A. P. Saunders Dr. H. Wroblowa 

The previously reported power densities for the porous I^/l electrode 
were underestimated. The main sources of error were due to air-locks in the 
back of the electrode and leaks around the electrode periphery, which pre- 
vented higher currents to be obtained at a high faradaic efficiency. 

The leaks have been a major retarding factor in this work due to the 
incompatibility of most materials with iodine. Graphite is the only material 
found to date suitable for the cathode, and two sealing materials, E 7 epoxy 
resin and RTV 60 silicon rubber have been used. The electrode (Fig. 1) has 
been used successfully for 13 days at a mean current density of over 200 
ma/cm . Examination showed that only the epoxy resin was visibly attacked, 
but this did not affect the cathode performance . 

This cathode is made of FC50 porous graphite (Pure Carbon Co.) with 
micro pores predominately in the 0.3-2 range ('median pore 0 ~ 0.8), and 
is combined with a Mg anode in a fuel cell. The oxidant is supplied to the 
porous cathode under pressure, and the current -potential characteristics are 
measured using model 610a and 150b Keithley electrometers and model #680 
Moseley recorders . 

Results 

Flow-rate measurements were carried out using an 0.2 cm thick and 
1 cm diameter graphite membrane. The flow rate, v, was measured by a de- 
crease in height of the meniscus of (0.973 ml 2 in 2 .5 m KJ) solution in a 
glass tube from which the electrolyte was pushed through the electrode under 
a pressure (Ap) . The change of height of the meniscus (Ah) was measured 
using a cathetometer . The time (At) was measured using a stop watch. The 
accuracy of the measurement was ^2.5%. Faradaic efficiency was calculated 

* Maximum current densities in excess of lA/cm^ have been drawn for 
periods of hours . 
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by comparing the observed flow rate with the simultaneously measured 
maximum constant current obtainable at the given Ap over a period of time. 
The efficiency was 100 ± 3% (i.e . ~100% within the experimental error) . 

The current -over potential relationship was recorded in steady state 
as a function of applied pressure. At each Ap value, there is only one 
point on the I curve corresponding to the condition 

I » 2F v C 
max o 


where I is the current and C is the concentration of I 9 in the supplied 
electrolyte, v is the flow -rate. At lower currents a part of non -reduced I 
passes into the solution at a rate measured by the value of (I max - I) . 

Discussion 

The current -potential relationship is shown .in Fig . 2 in form of 

plots of reduced current, , versus the cathode overpotential. The 

l max 

points seem to lie on the same curve for : > 0.9. This result is 

i max 

not predicted by theoretical analysis which would indicate an increase with 
the pressure gradient of the overpotential T| corresponding to the given 

4 value. However, for the range of Vp values studied, this increase 

l max . 

may be within the- limits of the experimental error. The numerical results 
of the theoretical analysis are being completed. 

On the basis of the results obtained, it may be anticipated that upon 
further increase of the pressure gradient, the increase of the cathode over- 
voltage corresponding to «95% 1^ consumption will not be very high. 

The power density of the Mg-I^ fuel cell obtainable in practical 
arrangement can now be projected. The power density, P, is given by 

P “ CV I2 - V Mg +1 ' I 2- 1R]i 

where V, • (at higher current densities) is approximately given by 
Mg 


V Mg = -1.45 + 1.67 i 

1 W. N. Carson, W. H. Fischer, E.'G. Siwek, Electrochem. Technology, 
5 , 423, 1967. 
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where V T = 0 . 535 v, TU (at 95% efficiency) - 0.35V, R ?=* 0.45 0 
L 2 l 2 

(specific resistivity of the solution is 4.5 Cl cm, distance between elec- 
trodes -0.1 cm) . 

Then 

P (95% eff) ~ [ 1,6 “ 2,12 1 - 1 1 

2 

and the maximum power density is obtained at i« 0.38 A/cm . The 

pressure gradient corresponding to this value is easily found from equation 

1 combined with the condition i = 0.95 i = 0.95 . 2F v C Vp, where v 

max o 

is the flow rate at Vp = 1 w2.2x 10 -6 l/sec, C = 1 M, is the initial 

cm o 

concentration of and Vp is the pressure gradient. 

Then 

VP . *»0.9 atm/cm 
opt 

The performance of the cell is limited mainly by polarization of the metal 
anode, which accounts for the major part of the power loss, and is the 
limiting factor of the cell operation. On the other hand, the power losses 
due to the cathode operating at about 0.95% efficiency of iodine consumption 
are small. 

In absolute values , the maximum power density of the operating cell 

2 

can be projected at P = 0.3 watts/cm . Further improvement would involve 
1) decrease of the anode polarization, 2) decrease of the IR drop in the 
electrolyte . 
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GAS BLEED 



Figure 1 
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Figure 2 


A3-11 




